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ABSTRACT 



Aims. To investigate the relationship between the X-ray variability amplitude and X-ray luminosity for a sample of 14 
bright Ultra- luminous X-ray sources (ULXs) with XMM- Newton /EPIC data, and compare it with the well established 
similar relationship for Active Galactic Nuclei (AGN). 

Methods. We computed the normalised excess variance in the 2-10 keV light curves of these objects and their 2-10 keV 
band intrinsic luminosity L2-10 kcv- We also determined model "variability-luminosity" relationships for AGN, under 
several assumptions regarding their power-spectral shape. We compared these model predictions at low luminosities 
with the ULX data. 

Results. The variability amplitude of the ULXs is significantly smaller than that expected from a simple extrapolation of 
the AGN "variability-luminosity" relationship at low luminosities. We also find evidence for an anti-correlation between 
the variability amplitude and L2--10 kcv for ULXs. The shape of this relationship is consistent with the AGN data but 
only if the ULXs data are shifted by four orders of magnitudes in luminosity. 

Conclusions. Most (but not all) of the ULXs could be 'scaled-down' version of AGN if we assume that: i) their black 
hole mass and accretion rate are of the order of ~ (2.5 — 30) x 1O 3 M0 and ~ 1 — 80% of the Eddington limit, and ii) 
their Power Spectral Density has a doubly broken power-law shape. This PDS shape and accretion rate is consistent 
with Galactic black hole systems operating in their so-called "low-hard" and "very-high" states. 

Key words. Black hole physics - X-ray: galaxies - X-ray: binaries 



1. Introduction 

Ultra-luminous X-ray sources (ULXs) are point-like sources 
with luminosities greater than 10 39 erg s _1 in the 0.3-10 
keV band. This high luminosity is greater than that ex- 
pected from stellar mass black holes (Mbh< 20 M©) ac- 
creting at the Eddington limit. Because they are usually 
located away from the nucleus of the galaxies they are 
unlikely to be associated with super massive black holes 
(SMBH, M B h> 10 5 M©), typically observed at the centre 
of active galactic nuclei (AGN). Their high luminosities 
can be explained if we assume that they host a black hole 
(BH) wit h an "intermediate ma s s" , ar ound 100-10000 M© 
(IMBHs, IColbert k MushotzkvL Il999h . However, the true 
nature of these sources is stil l unclear and oth er mecha- 
nisms as anisotropic emission (King et all [200l or accre- 
tion onto the BH in excess of the expected Eddin gton limit 
dBegelmanl . l2002j) could also explain them (see IRobertsl . 
l2007l fora recent review). The question of what powers 
ULXs will be conclusively answered by a direct mass mea- 
surement based on the determination of the binary orbit. 
However, due to their extragalactic nature, the study of the 
ULX counterparts in other bands has been difficult. 

In the meantime, both spectral and timing methods 
have been used over the last few years in an attempt 
to constrain the mass of the compact object in ULXs 



dMiller et al.l . l2004lh A common spectral method uses the 
temperature and luminosity of the accretion disc emission 
to determine the BH mass (assuming the standard Shakura- 
Sunyaev models). As for t he tim ing methods, one can either 
use the lM c Hardv etail (|2006l ) and iKording et all (|2007() 
scaling relationships of the characteristic time scales in 
Galactic black hole binaries (GBHs) and SMBHs with BH 
mass and accretion rate, or t he timing - spectral scalin 
for the QP Os in GBHs (e.g. IShaposhnikov k Titarchn " 
I2007ll2009t ). The results from the application of this method 
to ULXs have been non-conclusive. Som e studies suggest 
that ULXs are stellar mass BHs (e.g . iGladstone et al ' 



iGlax 

20091: IRobertsl . 120071 I Zezas et all 120071: iDewangan et al 
20061) while others imply that they host IMBHs (e.. 
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Casella et all 120081: iMiller et all . l2004jh 

In this paper we investigate the relationship between 
the X-ray variability amplitude and luminosity for a sample 
of 14 bright ULX sources using XMM-Newton /EPIC data. 
Our first aim is to measure the 2-TO keV normalised excess 
variance, cr^ xs , for these objects and investigate whether 
it correlates with the source luminosity. The normalised 
excess variance is a simple-to-calculate quantity that mea- 
sures the intrinsic variability amplitude of a source. It can 
be a useful complementary tool to the full-blown power- 
spectrum density (PSD) analysis, and it has the advantage 
that it can be applied to a larger number of objects as it 
does not require high-quality data (i.e. long, high signal-to- 
noise light curves). 
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Table 1. The sample and observational details. 



Name 


R.A. 


Dec 


Dist. 


ObsID 


Mode* 


Filter 


Back. Radius 


Seg.(T net 




(J2000) 


(J2000) 


(Mpc) 








(pixels) 


(ksec) 


NGC55ULX 


00 15 28.9 


-39 13 19.1 


1.9 


028740201 


FW 


Thinl 


534 


1(30) 


NGC253PSX-2 


00 47 32.9 


-25 17 50.3 


3.2 


152020101 


FW 


Thinl 


500 


2(40,35) 










125960101 


FW 


Medium 


500 


1(39) 


NGC1313X-1 


03 18 20.0 


-66 29 11.0 


1.0 


106860101 


FW 


Medium 


930 


1(31) 










405090101 


FW 


Medium 


677 


3(40,38) 


NGC1313X-2 


03 18 22.3 


-66 36 03.8 


4.0 


106860101 


FW 


Medium 


500 


1(31) 










405090101 


FW 


Medium 


500 


2(40,38) 


NGC2403X-1 


07 36 25.6 


+65 35 40.0 


3.5 


164560901 


FW 


Medium 


556 


2(40,30) 


HoIIX-1 


08 19 29.0 


+70 42 19.3 


3.3 


200470101 


FW 


Medium 


998 


1(37) 


M81X-6 


09 55 32.9 


+69 00 33.3 


3.7 


111800101 


SW 


Medium 


500 


2(40,37) 


M82X-1 


09 55 50.2 


+69 40 47.0 


4.0 


206080101 


FW 


Medium 


500 


2(40,37) 


HoIXX-1 


09 57 53.2 


+69 03 48.3 


3.3 


200980101 


LW 


Thinl 


737 


2(40,37) 


NGC3628X-1 


11 20 15.8 


+13 35 13.6 


12.0 


110980101 


EFW 


Thinl 


500 


1(38) 


NGC4559X-1 


12 35 51.7 


+27 56 04.1 


8.9 


152170501 


FW 


Medium 


1038 


1(37) 


NGC4945X-2 


13 05 33.3 


-49 27 36.3 


1.0 


204870101 


FW 


Medium 


500 


1(40) 


NGC5204X-1 


13 29 38.6 


+58 25 05.7 


5.3 


405690201 


FW 


Medium 


831 


1(37) 


NGC5408X-1 


14 03 19.6 


-41 22 59.6 


4.9 


500750101 


FW 


Thinl 


500 


1(40) 










302900101 


FW 


Thinl 


500 


2(40,40) 



POX52 12 02 56.9 -20 56 03.3 96.1 302420101 FW Medium 500 2(40,40) 



* 'FW, 'SW, 'LW and 'EFW stand for "Full", "Small", "Large" and "Extended Full Window" mode of the PN detector, respectively. 



It is well establish ed that g^yg i s anti - correlated with 
luminosity in AGN (jNandra et all Il997t iLeighlvl . Il999t 
iTurner et all Il999t) . More over, the excess variance anti- 
correlates with BH mass dLu fc Yul. 1200 It iBian fe Zhao| . 
2003 | :IO'Neill et all . 120051: iMiniutti et all 120091: IZhou et al 



20101) . Furthermore. iPapadakis! (|2004l) showed that the 



"variability-mass" relationship is probably the physically 
fundamental relationship rather than the "variability- 
luminosity" relationship in these objects. Our second aim 
is to compare the "variability-luminosity" relationship for 
ULXs with that of AGN with known BH mass. Following 
IPapadakis! (12001 . if we assume a universal PSD shape 
for AGN, which scales appropriately with BH mass and 
accretion rate, we can then make predictions on the ex- 
pected AGN "variability-luminosity" relationship at low 
luminosities. We want to investigate whether the ULX 
"variability-luminosity" data are consistent with various 
model "variability-luminosity" relationships for AGN, and 
if yes, what are the implications for the ULX PSD shape, 
BH mass and accretion rate. 

In Section 2 we present the sample selection. In Section 3 
we describe the data reduction, and in Section 4 we discuss 
the data analysis and present our results. We present a 
short discussion of their implications and our conclusions 
in Sections 5 and 6, respectively. 



2. The sample 

We considered all bright ULXs reported in the lite r ature , 



and in particular the object s studi ed bvlHei 



and IGladstone end] (|2009t) . The iHeil et all (|2009t ) sam 



et all (120091 ) 



pie includes all bright ULXs which have been observed 
with XMM-Newton for more than 25 ksec and their 0.2- 
10 keV flux is greate r than 5 x 10 _13 erg cm _2 s _1 . The 
IGladstone et al.l ( 20091 ) sample includes all ULXs observed 



with EPIC /XMM-Newton with more than 10000 net counts 
in the 0.3-12 keV EPIC band. 

One of our main aims is to compare the 2-10 keV vari- 
ability amplitude of ULXs with the variability amplitud e 
of the nearby AGN studied mainly bv lO'Neill etail (|2005l ). 
The length, T, and the bin size, At, of a light curve de- 
termines the lower and higher frequency sampled, since 
^min = 1/T and Vnmx = 1/(2 At) Hz. The excess variance 
of the light curve depends on the intrinsic power-spectrum 
and also on the minimum and maximum frequencies (see 
Section 4.5). For that reas on we used th e same length for 
the light curves as those of I O'Neill et al.l (|2005l) . Thus, we 
considered light curve segments with a length of 30 — 40 
ksec. Regarding At, due to the low count rate of all ob- 
jects in the sample, we used bins of size 1000 s in order to 
increase the signal-to-noise of their light curves. 

Consequently , we ch ose from the lHeil et al.l (|2009l ) and 
IGladstone et al.l (|2009() samples those sources which were 
observed by XMM-Newton with a net exposure, T net , larger 
than 30 ksec. For this reason we did not consider the XMM- 
Newton data of M33 X-8, IC342 X-l, NGC4395 X-l, and 
M 83 ULX. We did not consider NGC 4395 X-l either be- 
cause it was located on a gap of the PN detector during its 
XMM-Newton observation with T nct > 30 ksec. 

Our hnal sample comprises of 14 ULXs. Table [1] lists 
their coordinates, distance, and the XMM-Newton obser- 
vation details. Coordinates and distances were taken from 
the NASA/NEE0 database. Distances correspond to the 
average redshift-independent estimate for each object. For 
four sources (namely NGC 253 PSX-2, NGC 1313 X-l, 
NGC 1313 X-2, and NGC 5408 X-l) wc were able to re- 
trieve from the archive two observations with T not larger 
than 30 ksec. Note that the starburst galaxy M82 con- 
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tains two ULXs, namely X41.4+60 and X42.3+59, which 
are unresolved by XMM-Newton, and they both contribute 
to the M82 X-l light curve. During the 2004 April ob- 
servation that we considered in this work, approximately 
84% of the observed c ount rate originates from X41.4+60 
tiFeng fc KaaretL 12001 . 

In order to extend the lO'Neill et al.l (|2005f ) sample to 
include AGN with low BH masses, we also considered the 
XMM-Newton observation of POX 52, which h osts an AGN 
with a low BH mass (M BH = 1.6 x 10 5 M Q ; iBarth et all 
l2004t) . The coordinates, distance, as well as the XMM- 
Newton observation details for this source are also listed 
in Table Q] (the distance in this case corresponds to the 
"luminosity distance" estimate of NASA/NED). 



3. Data reduction 

Data were retrieved from the XMM-Newton public data 
archival. We used the XMM-Newton Science Analysis 
System SA£0 software version 9.0.0 and followed stan- 
dard procedures to extract science products from the 
Observation Data Files (ODFs). 

We used data from the EPIC-pn camera only due to 
its superior statistical quality. Source counts in each case 
were accumulated from a circular region of radius 400 pix- 
els, centred on the source's RA and Dec. In the case of 
NGC253 PSX-2 (ObsID 152020101) we used a radius of 
300 pixels to avoid contamination from a nearby source and 
the detector gap. Background data were extracted from a 
source free circular region on the same CCD chip than the 
source (background region radii are listed in Column 8 of 
Table [I]). We selected only single and double pixel events 
(i.e. patterns of 0-4). Bad pixels and events too close to the 
edges of the CCD chips were rejected using "FLAG=0." 
Given the observed count rate, photon pile-up is negligible 
for the PN detector in all cases. 

Source and background light curves in the 2-10 keV 
band were extracted using evselect task on SAS with 
a 1000-sec bin. They were screened for high background 
(usually at the end and/or the beginning of the individual 
observations) and flaring activity. After rejection of the re- 
spective time intervals, the total useful observation time for 
each observation is usually less than the original PN expo- 
sure time. We chose to study only those light curves with 
at least one "clean" segment longer than 30 ksec. 

As an example, in the top panel of Fig. [T] we show the 
light curve of NGC 5408 X-l, which is typical of the light 
curves of all sources in the sample. Filled circles and stars 
indicate the 2-10 keV (background subtracted) source and 
background light curves, respectively. The background light 
curve corresponds to the full PN exposure length while 
the source light curve is plotted only for those parts of 
the observation when the background activity was "low". 
In general, as background "loud" we identified the obser- 
vation parts where: i) the background light curve showed 
"flare"-like events and/or prominent decreasing/increasing 
trends (usually at the start/end of an observation), and ii) 
the "net" source count rate was less than twice the back- 
ground count rate. The brackets on top of the NGC 5408 
X-l light curves indicate the light curve segments we chose 
for this observation. Clearly, the background light curve is 
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Fig. 1. Light curves of NGC 5408 Xl/ObsID 302900101 
(top panel) and POX 52 (bottom panel). Stars indicate the 
background light curves and filled dots indicate the back- 
ground subtracted light curves, plotted only for the ob- 
servation period for which the background activity is "low" 
(see text for details). The brackets on top of the light curves 
indicate the segments that were used to estimate the excess 
variance in each case. 



2 http:/ /xmm.esac.esa.int/xsa/index.shtml 

3 http://xmm.esac.esa.int/sas/ 



stable, and of much less intensity than the "net" light curve 
count rate. This was the case for almost all of the light 
curve segments we used in this work. We also extracted the 
EPIC-pn spectra, after the rejection of the time intervals 
affected by high background, using single and double events 
(PATTERN <=4). Response and auxiliary matrices were 
created with SAS tools RMFGEN and ARFGEN, respectively. 

The POX 52 XMM-Newton data were reduced in the 
same way. Its 2-10 keV (background subtracted) source 
and background light curves are also plotted in Fig. [1] (bot- 
tom panel) . The observation is affected by high background 
flaring activity during the first ~ 20 ksec, and after ~ 40 
ksec since the start of the observation, which lasted for al- 
most 10 ksec (note that this observation shows the worse 
background "flaring" activity among all the light curves we 
studied in this work). 



4. Data analysis and results 

4.1. The variability amplitude estimation 

As a measure of the intrinsic variability amplitude of the 
light curves we computed their normalised excess variance, 
°nxs- ^ s sc L uare root is a measure of the average variability 
amplitude of a source as a fraction of the light curve mean. 
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where x, cr err , and N are the count rate, its error, and the 
number of points in the light curve, respectively and S 2 is 
the variance of the light curve, i.e.: 

1 N 

i=l 

Following lO'Neill et al.l (l2005h . we computed a^xs f° r eacn 
continuous light curve segment with 30 < T nct < 40 ksec. 
For light curves longer than 40 ksec we considered only the 
first 40 ksec. If there were more than one segments of dura- 
tion T ne t larger than 30-40 ksec we computed er NXS f° r eacn 
one of them. The number of light curve segments in each 
observation, and their T not , are listed in Table [1] (Column 
9). A few "missing" points within each segment, due to 
the presence of background flaring activity, appear in one 
of the light curve segments of NGC4559 X-l, NGC4945 
X-2, NGC253 PSX-2 (ObsID 152020101), and HoII X-l. 
Missing points are typically less than 10-15% of the total 
number of points. The first segment of POX 52 light curve 
shows the largest number of "missing" points (20% of the 
total). The presence of missing points in these segments 
should increase the uncertainty of the resulting crj^xs esti- 
mates. 

Our cr NXS estimates, for each light curve segment, are 
listed in Table [2] (Column 3). The numbers in parenthesis 
in Table [5] indicate the weighted mean <7 NX s an< ^ ^ s error 
in the case we had more than one excess variance estimate 
for the same source. 

For two sources the ct^xs measurement was negative. 
In these cases, we estimated the 90% upper limits of the 
intrinsic o" xxg values using the 90% upper limit s on the 
source variance, as listed bv lVaughan et al.l (|2003|) in their 
Table 1. In order to constrain as much as possible the 
"variability-luminosity" ULXs correlation (see Section 4.3 
below), we assumed a PSD slo pe of —1, and the upper limit 
from the lVauehan et aTl (|2003h simulations with the longest 
light curves (any other choice would result to an even larger 
90% limit). In order to take into account the uncertainty on 
(T NXS due to the experimental Poisson fluctuations as well, 
we added to these limits the value of 1.282 err(crf5[ XS ). Our 
final estimates of the 90% confidence limits for these two 
sources are listed in a par enthesis in Tabl e [2| 

Regarding M82 X-l, iFeng fc Kaaretl (|2007ft estimate 
that X41.4+60 contributes more than 80% of the observed 
count rate during the 2004 April XMM-Newton observa- 
tion of M82. X42.3+59 is high ly variable source, bu t on 
time scales of years (see Fig. 5 in lFeng fc Kaaretl I2007T ). On 
shorter time scales, the same authors show that the PSD 
of X41.4+60 has a significantly higher amplitude than the 
PSD of X42.3+59. 

4 The error accounts only for the uncertainty due to the 
Poisson noise and not due to the stochastic character of the 
intrinsic variability process. 



Table 2. The excess variance (<r NXS ) and 2-10 keV intrinsic 
luminosity (L2-10 kov) in logarithmic scale. 



Name 


ObsID/Seg. 


°NXS 


log(Lx) 






(xlO -3 ) 




NGC55ULX 


028740201/1 


147 ± 2.1 


38.7 


NGC253PSX-2 


152020101/1 


5.7 ±2.4 


39.3 




/2 


-1.8 ±2.6 






125960101/1 


10.0 ±2.0 


39.4 




(6.5 ± 1.3) 


(39.4) 


NGC1313X-1 


106860101/1 


2.4 ± 1.6 


39.6 




405090101/1 


4.0 ± 1.5 


39.6 




/2 


-0.3 ±1.5 








(2.0 ±0.9) 


(39.6) 


NGC1313X-2 


106860101/1 


-7.1 ±6.8 


39.1 




405090101/1 


22.9 ± 1.4 


39.7 




/2 


11.3 ± 1.3 






(16.2 ±0.9) 


(39.5) 


NGC2403X-1 


164560901/1 


-2.3 ±3.5 


39.2 




/2 


-2.8 ±4.5 






(< 8.3) 




HoIIX-1 


200470101/1 


1.3 ± 1.0 


39.6 


M81X-6 


111800101/1 


0.9 ± 2.2 


39.5 




/2 


8.2 ± 2.1 






(4.9 ± 1.5) 




M82X-1 


206080101/1 


1.0 ±0.3 


40.4 




/2 


0.8 ±0.3 






(0.9 ±0.2) 




HoIXX-1 


200980101/1 


0.1 ±0.6 


39.8 




/2 


1.2 ± 0.6 






(0.6 ±0.4) 




NGC3628X-1 


110980101/1 


0.6 ± 6.2 


40.0 


NGC4559X-1 


152170501/1 


13.0 ±6.6 


39.6 


NGC4945X-2 


204870101/1 


-10.5 ±9.5 


39.0 






(< 27.0) 




NGC5204X-1 


405690201/1 


1.6 ± 2.5 


39.6 


NGC5408X-1 


500750101/1 


16.0 ±5.1 


39.3 




302900101/1 


2.4 ±5.0 


39.5 




/2 


12.4 ±5.1 






(10.2 ±3.0) 


(39.4) 


POX52 


302420101/1 


78 ±17 


40.7 




/2 


105 ± 15 








(93 ±11) 





4.2. The hard band X-ray luminosity estimation 

To estimate the X-ray luminosity for each source we fit- 
ted their spectra with an absorbed power-law model in 
the 2-10 keV band. For the Galactic absorption, we fixed 
the Nh values at th e value s derived from the HI maps of 
(jDickev fc Lockmanl I1990D . The spectral fitting was per- 
formed using XSPEC version 12.5.1. Using the best-fit re- 
sults we estimated the source flux in the 2-10 keV band, and 
hence the source luminosity, L2-10 kcVi adopting the dis- 
tance estimates listed in Tablc[TJ The unabsorbed L2-10 keV 
estimates are listed in Table [2] (Column 4)0- The values in 
parenthesis correspond to the mean L2_io keV estimates, in 
the case there were more than one spectrum for an object. 

J Note that we considered a single spectrum for each observa- 
tion (i.e. we accumulated all the data for the whole "background- 
quiet" period of each observation), irrespective of the number 
of segments that we used for the estimation of the normalised 
excess variance. 
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38.5 39.0 39.5 40.0 
Log[L x (2-10keV)] 



40.5 



Fig. 2. Normalised excess variance versus log(L 2 _io kov) 
for the ULXs. Arrows indicate the 90% confidence upper 
limits on ct^xs for the sources with negative excess vari- 
ance estimates. The point plotted with an empty square 
indicates the M82 X-l measurement when corrected for 
the contribution of X42.3+59 to the observed count rate 
(see text for details). The shift of M82 X-l is shown as a 
dotted line. The solid line indicates the best-fit to the data 
(excluding NGC55 ULX). 



In half of the cases, the best-fit x? e d values were larger 
than ~ 1.2. This is mainly due to the presence of additional 
complexity in the spectra that the simple power-law model 
cannot account for (e.g. the presen ce of "breaks" in the hig h 
energy spectra of these sources, iGladstone et all . 120091) . 
Nevertheless, the power law model describes adequately the 
broad shape of the source spectra in all cases, and the re- 
sulting best-fit flux measurements should be an accurate 
estimate of the source X-ray continuum flux. To investi- 
gate this issue further, w e used the 2-10 keV band best-fit 
results of IStobbart et all (|2006fl to estimate the 2-10 keV 
luminosity for the nine sources in common. We found that 
Lours = Literature m ah cases except for NGC 55 ULX, Ho II 
X-l, and Ho IX X-l sources where L ours /Liiteraturc = 1-1, 
Lours/Liiterature = 0.8, and L ours /Literature = 0.8, respec- 
tively. We are thus confident that o ur luminosity est i mates 
are reliable. Regarding M82 X-l, iFeng fc Kaaretl (|2007t) 
estimate a 2-10 keV luminosity of 1.7 x 10 40 ergs s~ 4 for 
X41.4+60 during the April 2004 XMM-Newton observation 
of the source. This is smaller than our estimate of 2.5 x 10 40 
ergs s _1 , but this is expected given the presence of the other 
ULX, which also contributes to the flux we measure from 
this source. 



4.3. The "variability-luminosity" relation of ULXs 

Fig. [2] shows ct^jxs as function of log(L2_io kev) for the 
ULXs in the sample. The arrows indicate the 90% con- 
fidence upper limits on the (intrinsic) excess variance of 



the two sources with negative cr^xs estimates. The X-ray 
variability amplitude appears to decrease with increasing 
X-ray luminosity. To a large extent, this trend is driven by 
the NGC 55 ULX data. However, this source shows "dip- 
ping" episodes in its variability , whic h enhance its variabil- 
ity amplitude ijStobbart et all |2004|) . Similar events have 
not been observed in other ULXs. 

To investigate the significance of the apparent 
"variability-luminosity" relation in Fig. [5] we fitted the 
[log(er^ xs ), log(L 2 -io kev)] data with a straight line of the 
form log(o-Nxs) = a+bxlog(L 2 _io kev)- Given the presence 
of upper limits in two objects, we used the Buckley- James 
regressio n method as imple mented in the software package 
ASURV (llsobe et al.l . ll986h . Since the NGC 55 ULX vari- 
ability properties may be somewhat "anomalous" among 
ULXs, we excluded this source from the fit. 

The best-fit slope value is a = —1.0 ± 0.4 and is sig- 
nificantly different from zero at the 2.5<r level. The best- 
fit is indicated by the solid line in Fig. O The point plot- 
ted with an empty square in Fig. [2] indicates the M 82 X- 
1 measurements, when "corrected" for the contribution of 
X42.3+59 to the observed count rate . We have adopted the 
2-10 keV luminosity measurement of lFeng fc Kaaretl (|2007f ) 
for X41.4+60. Furthermore, we have increased our excess 
variance measurement by a factor of 1.5. This is based 
on the fact that the "area A" and "area A+B" PSDs of 
IFeng fc Kaaretl (|2007h show a flat PSD at low frequencies, 
whose normalisation is higher by <~ 1.5 in the case of "area 
A" PSD. In the following figures we indicate only the "cor- 
rected" data for M82 X-l, albeit with a different symbol 
than the rest of the ULXs data. 



4.4. Comparison with AGN 

The top panel in Fig. [3J shows the "variability-luminosity" 
plot for the ULXs (filled squares in all panels, except for the 
M 82 X-l data which are shown with an open square; ar- 
rows indica te 90% confidence up per limits) and AGN, using 
the data of lO'Neill et all ()2005l ) (open triangles). We also 

data for the 4 



added to this plot the [( Tmyc; , log(L2-inkeV J I u i 
AGN with IMBHs from IMiniutti et all (|2009h (black aster- 
isks). The open circles in all panels indicate the POX 52 
data. The variability amplitude of this source is compara- 
ble to the amplitude of th e two lowest luminosity objects 
in the lO'Neill et al.l (|2005l ) sample a nd with the amplitude 
of the lowest luminosity object in the IMiniutti et aiT(|2009|) 
sample. In fact, the addition of the POX 52 and the four 
IMBHs data in the plot strengthens the possibility that the 
AGN "variability-luminosity" relation may flatten at lumi- 
nosities lower than ~ 10 42 ergs s _1 . 

The solid line in the top panel of Fig. [3J indicates the 
best-fit line for the ULXs data while the dashed line indi- 
cates the same line shifted by A log(L2_iokcv) = 4. The 
shifted line describes rather well the AGN "variability- 
luminosity" relation. In fact, we used th e "ord inary least 
squares bisector" method of llsobe et al.l (|1990l ) to fit the 
AGN data (in the log-log space) with a straight line. The 
best-fit slope was —1.22 ± 0.16 which is consistent with 
the best-fit slope for ULXs (&agn — ^ulxs = 0.2 ± 0.4). 
Therefore, the variability amplitude may indeed decrease 
with increasing luminosity in a similar way for AGN and 
ULXs. 

Despite this similarity, Fig. [3J also indicates that the 
ULX data are not consistent with the AGN data. For a 
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given luminosity, even if the AGN "variability - luminosity" 
relation flattens below ~ 10 42 ergs s _1 , the ULX variability 
amplitude is at least 10 times lower than expected when 
we extrapolate the AGN "variability-luminosity" relation 
to lower luminosities. Only NGC 55 ULX appears to be 
consistent with the AGN data. 



4.5. Determination of model 
relations 



variability-luminosity" 



To compare in a quantitative way the ULXs and AGN 
"variability-luminosity" relationship it is necessary to de- 
rive the "excess variance - luminosity" relationship for AGN 
at low luminosities, and then compare it with the observed 
relationship for the ULXs. In this way we will be able to 
investigate if and for which physical parameters (i.e. MbHi 
and accretion rate in units of the Eddington limit, m E dd) 
the ULXs data will be in agreement with the "model" AGN 
"variability-luminosity" relationships . 

The bolometric luminosity emitted by an AGN is Lboi = 
m E dd^Edd = 1.3m E ddl0 38 (M B H/MQ ergs s" 1 ). If k bo i is 
the X-ray to Lboi conversion factor then, 



L 2 -iokoV = fc bo; 1.3m E ddlO (M B h/M ) ergs s" 



(4) 



(in all the calculati ons below, we adop ted the kb i — L E dd re- 
lationship given by iLusso et all . [20091) . On the other hand, 
the observed excess variance, estimated from a light curve 
of length T and bin size At, is an approximate measure of 
the integral: 



'NXS 



P(i/)di/, 



(5) 



where P(^) is the intrinsic power spectrum normalised 
to the square of the light curve mean, i/ m i n = 1/T and 
i/ max = 1/(2 At). Equation (5) results from Parseval's the- 
orem for Fourier series, and from the fact that the mean 
value of the periodogram (i.e. the square of the discrete 
Fourier transform of the light curve) is (approximately) 
equal to P(^) (e.g. s ee discussion given in Section 2.2 of 
IVauehan et all I2003T ). 

The excess variance can be associated with the BH mass 
and accretion rate if one assumes a PSD shape and certain 
scaling relations between the PSD characteristic frequencies 
with Mbh and m E dd ■ Below we present model "variability - 
luminosity" relations for AGN, assuming two rather simple 
scenarios for their PSD shapes, which are based on recent 
power spectral studies of AGN and GBHs. 

Case A : Analyses based on high quality RXTE and 
XMM-Newton light curves have shown that the AGN 
PSDs can be approximated by a broken power law, with 
~P(v) = A^/v^Y 1 and V(v) = A(v/v hl )~ 2 at frequen- 
cies below and above a charact e ristic frequency break , 
yy,r fe.g. iMarkowitz etahl 120031: lM c Hardv et"aU 12004ft . 
I M c Hardy et al.l ( 20061 ) have demonstrated that z-v depends 
on both A/bh and m E dd as: 



Vbr = 0.003m Ed d(M B H/10 6 M Q )" 



((») 



In this case (case A model hereafter), it is straightforward 
to show that: 



'NXS 



CiH^)-^" 
C i ln (» 



(if Z/b r < I/min) 

(if Z/ m in < V hr < V-n 

(if V hr > U max ) 




41 42 
Log[U(2-10 keV)] 

Fig. 3. Normalised excess variance versus log(L 2 _iokcv) 
for ULXs (black squares, and open square for the M82 
X- l data) and AGN (o pen triangles). IMBHs reported 
by iMiniutti et all ([20091) are included as black stars and 
POX 52 (reported here) as an open circle. Top panel: the 
solid line indicates the best fit to the ULX data, and the 
dashed line indicates the same line shifted by +4 along 
the x-axis. Middle panel: The lines indicate the case A 
model "variability-luminosity" relation for m E dd = 0.03 
(dotted line), 0.1 (dashed line) and 0.3 (continuous line). 
The short solid lines between the model curves indicate 
black hole masses of 5 x 10 6 , 5 x 10 7 , and 5 x 10 8 M Q , from 
left to right. Bottom panel: The lines indicate the case B 
model "variability-luminosity" relation, for the same accre- 
tion rates. As above, the short solid lines within the model 
curves indicate black hole masses of 10 3 , 10 4 , 5 x 10 5 , 5 x 
10 6 , 5 x 10 7 , and 5 x 10 8 M , from left to right. 



where Ci = A^br- Following iPapadakisl (|2004f ) we assumed 
that C\ = 0.02. The solid line in Fig. 0] indicates the case 
A PSD model, when v ht = 10~ 3 and A = 20 Hz~ x (values 
chosen arbitrarily). The vertical, dotted lines indicate the 
v-min and i/ max frequencies. The diagonal dotted lines indi- 
cate the area below the case A PSD curve, and between v m i n 
and zAnax values, which according to equation (5) should be 
(approximately) equal to cr^xs- 

Case B : A second frequency break, Vbr,2i below which 
the PSD is roughly flat (i.e. P(z^) oc v°), has also been ob- 
served in at least one AGN (i.e. Ark 564, iPapadakis et all 
[200l lM c Hardv etHl [2001 and in GBHs in "low/hard" 
and "very high" states (see discussion in Section 5.1 below). 
The PSDs in the latter case are quite c omplex, usually de- 
scrib ed by a series of Lorentzians (e.g. iPottschmidt et all 
120031 ). However, the entire spectral shape roughly resem- 
bles a (doubly) broken power law of the form: P(is) = 
(7)A(v/v hl )~ 2 , for v > v hr , P(i/) = A^/v^)- 1 , for zv, 2 < 
v < nubr and P(f) = ^(^br^/^br) -1 =constant, for 
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Frequency (Hz) 



Fig. 4. The assumed PSD shape in case A (continuous line) 
and case B (dashed line) scenarios discussed in Section 1431 
Vbr and Vbr,2 are the characteristic "frequency breaks" , at 
which the PSD slope changes, while v min and v max are 
the minimum and maximum frequencies sampled by the 
light curve segments used in this work. The regions filled 
with diagonal-dotted lines (vertical-continuous lines) indi- 
cate the area below the PSD shape in case A {case B) sce- 
nario, which should be approximately equal to the excess 
variance of the light curves. This area (and hence er^xs 
as well) depends on the location of the break frequencies 
(determined by black hole mass and accretion rate) with 
respect to i/ max and v m i n (which are fixed according to the 
values T and At of the light curves used in this work). 

v < fbr,2- hi this case (case B model hereafter), the ex- 
cess variance of the light curves should be: 

(if J/br,2,14>r < fmin) 
(if fmin < t'br < "max and l^ b r,2 < "min) 
(if I/ br > fmax and ^ br ,2 < U min ) 

Ci[ln(^-) +2 - ^ - -a*-], 

(if Umin < V hl - 2 < Vbr < "max) 



a NXS — ' 



(8) 



C 1 n n (^) + 1-^1, 

L V "br,2 ' ^br,2 J ' 

(if "min < "br,2 < "max and t/ br > Z^nax) 
Cl(f max — "min)/"br,2, 

(if f br ,2," b r > fmax) 

Th e ratio Vhrlvhr.-i is usu ally ~ 10 in GBHs. We adopted 
the lAxelsson et al. | (|2006l ) relation between the two break 



(9) 



frequencies in Cyg X-l q 

^br,2 = 0.15^ r 2 . 

6 Strictly speaking, equation (9) is valid for Cyg X-l in its 
"low/hard" state. We verified that our results do not depend 
strongly on the adopted relationship between Z4, r ,2 and i-v, as 
long as Vhi,2 ~ 0.2 — 0.05^- 



The dashed line in Fig. H] indicates the case B PSD model, 
assuming a second break frequency which is 10 times 
smaller than j/ br . The vertical lines between v m i n and i/ max 
indicate the area below the the case B PSD curve, which 
should be approximately equal to cf^xs- We can now use 
the equations above to construct AGN model "variability- 
luminosity" relations as follows. 



4.6. The AGN model "variability-luminosity" relations 

We considered BH mass values in the range between 10 3 — 
10 9 M©, and three accretion rate values, namely 0.03, 0.1 
and 0.3 of the Eddington limit. For any given Mbh and 
"T-Edd values we used equation (4) to compute the 2-10 
keV luminosity of the source, and equations (6) and (9) to 
compute z^ br and f br .2, respectively. We then used equations 
(7) and (8) to compute the model cr^xs va hies in case A and 
case B scenarios, respectively. The minimum and maximum 
sampled frequencies in our case are z/ max = 1/(2 x 1000) 
Hz and i/ min = 1/T mca n Hz, where T moan = 37 ksec, i.e. 
the average length of all the segments listed in column (9) 
of Table 1. The resulting model "variability-luminosity" 
relations are plotted in Fig. [3] (middle and bottom panels). 

The dotted, dashed and solid line in the middle panel of 
Fig. [3] indicate the expected case A "variability-luminosity" 
relations for riiEdd =0.03, 0.1, and 0.3, respectively. The 
agreement between these lines and the AGN data is reason- 
ably good, indicating that both the shape and the scatter 
in the observed "variability-luminosity" relation for AGN 
can be explained if the nearby, bright Type-1 Seyfcrts ac- 
crete at ~ 3 — 30% of the Eddington limit. The flatten- 
ing of the relation at low luminosities is due to the fact 
that, if equation (6) is valid, then i\, x > u maK for objects 
with BH mass (luminosity) less than ~ 2 - 6 x 10 5 M Q 
[(0.6 — 4.0) x 10 41 er<7 s -1 ]. Consequently, the excess vari- 
ance should remain constant [sec bottom relationship in the 
set of equations (7)] for all objects with smaller Mbh (and 
therefore source luminosity). 

The dotted, dashed, and solid line in the bottom panel 
of Fig. [3J indicate the case B "variability-luminosity" pre- 
dictions (as before, we considered the values of riiEdd =0.03, 
0.1, and 0.3, respectively). If there is indeed a second PSD 
frequency break, and equations (6) and (9) are valid, then 
we expect fb r ,2 > "min for sources with BH mass (lumi- 
nosity) smaller than ~ 6 x 10 5 M Q (~4x 10 41 ergs s _1 ). 
Since the low-frequency flat part of the PSD does not con- 
tribute to the integral defined by equation (1) as much 
as the v~ x part does (see Fig. 0]), the excess variance 
is expected to decrease with decreasing luminosity when 
L 2 -iokcV < 4 x 10 41 erg s _1 . 



4.7. The comparison between ULXs and AGN revisited 

It is clear from the middle panel of Fig. [3] that, apart from 
NGC 55 ULX, the ULXs data are not consistent with the 
case A AGN model predictions. On the other hand, the bot- 
tom panel in Fig. [3J shows that most of the ULXs data are 
consistent with the case B model "variability-luminosity" 
relations. 

In Fig. [5] we plot again the <T^ xg — log(Zv2-iukev) data 
for the ULXs. The dashed lines in the same plot indicate 
the case B model "excess variance - luminosity" relations 
for AGN-like objects with a BH mass of (2.5 - 30) x 10 3 
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Fig. 5. Normalised excess variance versus log(L 2 _io kov) 
for the ULXs. The dotted-fillcd region indicates the area 
with BH mass of (2.5 — 30) x 10 3 M Q and accretion rates 
of mEdd=(0.01-0.8). Dashed lines indicate the case B model 
for AGN-like objects with a BH mass of (2.5 - 30) x 10 3 
M Q (from left to right along each line), and for accretion 
rates of m Ed d=[0.01, 0.05, 0.1, 0.2, 0.5, 0.8] (from top to 
bottom). 



Mq (from left to right along each line), and for various 
accretion rates. The dashed line on the top-left corner of 
the plot corresponds to mEdd = 0.01. Moving down along 
the lines, the accretion rate increases up to mEdd = 0.8. 
The fact that most of the ULXs data are located within 
the boundaries of the diagonally dashed region of the plot 
implies that these objects may operate like AGN, with a BH 
mass in the range (2.5 — 30) x 10 3 M Q , and an accretion 
rate in the range [0.01-0.8]. Regarding NGC 2403 X-l 
and NGC 4945 X-2 (the upper limits) , their low variability 
amplitude can be explained if they host a BH with a mass 
close to or even lower than 2.5 x 10 3 M Q . 

Furthermore, the same figure can also explain the ap- 
parent ULX "variability-luminosity" anti-correlation. If the 
higher-luminosity objects in the sample have systemati- 
cally larger BH mass and accretion rate than the lower- 
luminosity objects, then their variability amplitude should 
also be systematically smaller, hence the "smaller variabil- 
ity amplitude with increasing luminosity" trend we de- 
tected. 

Finally, it is clear from Fig. [3] that NGC 55 ULX is not 
consistent with the case B model predictions. One could 
assume that this source is more consistent with the case 
A model predictions, however, the high amplitude dipping 
episodes seen in the light curve of this source are not com- 
monly seen in AGN light curves. Furthermore, if that were 
the case, we would expect that ^b r > v max for this object, 
and its PSD to have a —1 sl ope in the freque ncy range 
- 10~ 4 - 10~ 3 Hz. However, iHeil et all (|2009D found an 
average PSD slope of —1.96 ± 0.04 in this frequency range. 



Therefore, the agreement of the NGC 55 ULX data with 
the case A "excess variance - luminosity" relation must be 
coincidental. 



5. Discussion 

We present the results of a variability analysis of a sam- 
ple of 14 bright ULXs using 19 observations with XMM- 
Newton /EPIC. We calculated their normalised excess vari- 
ance using light-curves of 40 ksec length. Our main aim was 
to compare their "variability-luminosity" relationship with 
the same relationship for AGN. Our main results can be 
summarised as follows: 

— The variability amplitude of ULXs is significantly 
smaller than that expected from a simple extrapola- 
tion of the AGN "variability-luminosity" relationship 
to lower luminosities. 

— We found evidence that the variability amplitude in 
ULXs decreases with increasing 2-10 keV source lumi- 
nosity. This "variability - luminosity" anti-correlation 
is similar (in slope) to what is observed in nearby Type 
1 Seyferts. 

We discuss in some detail some implications of our re- 
sults below. 



5.1. Are most ULXs AGN-like objects? 

The fact that ULXs show a significantly smaller variability 
amplitude (when compared to the amplitude expected from 
an extrapolation of the AGN "variability-luminosity" rela- 
tion to low luminosities) is consistent with the hypothesis 
that ULXs arc "scaled-down" version of the nearby AGN, 
but only if: a) there two break frequencies in their PSDs, b) 
Vbr scales with BH mass and accretion rate as in equation 
(6), c) ^br/^br,2 ~ 5 — 50, and d) 24,^2 (i.e. the frequency 
where the PSD shape changes from a slope of to ~ —1) 
is higher than 3 x 10 -5 Hz. 

This band-limited noise PSD shape is commonly 
observed in GBHs in their "low/hard" and "very 
high" state (VHS; see for exa mple Fig. 4e and 4c in 
iKlein-Wolt fc van der Klisl . 120081 ). In this case the variabil- 
ity is mostly limited to ~ 1 — 2 decades of temporal fre- 
quency. In most AGN though, the PSDs resemble that of 
GBHs in their "high state" : the 1/v part of the power 
spec trum extends over m any decades of frequency below 
u h[ (|M c Hardv et all . 120061) . The only exception is Ark 564, 
where a second PSD break (to a slope flatter than —1) is 
observed, and most of the variability is indeed limited in 
less than two decades of frequency. Interestingly, Ark 564 is 
the only high accretion rate AGN for which a good qual- 
ity PSD is available at the moment. lM c Hardv et all (|2007f) 
suggested that this object is the AGN analogue of GBHs in 
VHS. Perhaps then, just like Ark 564, ULXs are also similar 
to GBHs in VHS. 

Are there any indications that ULXs power spectra have 
a shape consistent with t he PSD model outlined above? 
Recently, IHeil et aD (|2009h presented the results from a de- 
tailed PSD analysis of archival ULX light curves. They de- 
tected "0 to -1" PSD breaks at frequencies higher than 
3 x 10~ 5 Hz in NGC 5408 X-l and M 82 X-l. In NGC 1313 
X-l they detected a "0 to -2" break, however the best-fit 
high-frequency slope has a large uncertainty (—2.35 ± 1). 
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The possibility of a -1 slope cannot be excluded, as it is 
just 1.35cr away from the best fit value (note also that "0 
to -2" PSD slope breaks are unusual in GBHs). 

On the other hand, the same authors found that the 
PSD of Ho IX X-l, NGC 1313 X-2, and NGC 55 ULX were 
best-fitted by a simple power-law model. The —2 PSD slope 
in NGC 55 ULX is not consistent either with the case A or 
case B scenario (see also Section |4~7]) . The best- fit slope for 
NGC 1313 X-2 is also steep (~ —2), but it is ba sed on just 
one p oint (see e.g. the top left panel in Fig. 1 of iHeil et all 
l2009t) . while the best-fit slope of Ho IX X-l is flat, close to 
~ —0.5. This may be indicative of a Vbv,2 break just below 
the lowest frequency of 10 -4 Hz that lHeil et al.l (l2009f) con- 
sidere d. We therefore believe that the results of IHeil et al.l 
(|2009f) do not disagree with the case B PSD shape we out- 
lined in the previous section. 



5.2. Constrains on the BH mass and accretion rate of ULXs 

If indeed most of the ULXs in the sample operate like AGN, 
and their PSDs are consistent with the case B scenario 
outlined in Section 4.5, then they must host a black hole 
with a mass of Mbh ~ 2500 — 3 x 10 4 M Q , which accretes 
at <~ 1 — 80% of the Eddington limit. Smaller BH masses 
with a higher accretion rate, or higher BH masses with a 
smaller accretion rate cannot be consistent with most of the 
ULX data shown in Fig. [3] To illustrate this point, we no- 
tice the solid line for 10 3 M Q in the bottom panel of Fig. [31 
which indicates the expected [of ms , log(L 2 _iokcv)] relation- 
ship in the case B model for objects with M B h = 10 3 M 
and 0.03 < niEdd < 0.3. Clearly, this line is not consistent 
with the ULX data. An increase of kbol by a factor of ~ 10 
would be required to shift this line along the x-axis by the 
same factor, and hence to be roughly consistent with the 
ULX data. This difference in kbol would imply a significant 
difference between the ULX and AGN spectral energy dis- 
tributions. Thus, if indeed most of the ULXs in our sample 
host a BH mass significantly smaller than (1 — 10) x 10 3 
M Q , then they are not exactly like AGN. 

IWu fc Gul (|2008f ) also addressed the questions "do ULXs 
operate like- AGN?", and, "if yes, what should their BH 
mass be"? They studied the relationship between T and 
X-ray luminosity in seven ULXs (four of them included in 
our sample). They found that it is similar to what is ob- 
served in GBHs and AGN, and that the ULX central BH 
mass should be ~ 10 4 M^ , in ag reement with our results. 
IStrohmaver fc Mushotzkvl (|2009f ) showed that the pattern 
of spectral and temporal correlations in NGC 5408 X-l is 
analogous to that seen in GBHs, and argued that the BH 
mass range for this system is from (2 — 9) x 10 3 M Q . The 
position of the NGC 5408 X-l data in Fig. [2] is consistent 
with an object with a BH mass of ~ 6 x 10 3 M Q which ac- 
cretes at ~ 4% of the Eddington limit, in agreement with 
the BH mass range of IStrohmaver fc Mushotzkvl (|2009h . As 
for the M 82 X-l data in the same figure, they suggested a 
BH mass of 1.8 x 10 4 M Q in the system. This is 2-3 times 
larger than the most likely BH mass range of "one to a 
few thousand solar masses" for this source, as estimated by 
iKaaret et all (|2009f ). This is not a significant discrepancy, 
because we would need many observations to estimate the 
average variability amplitude and X-ray luminosity of the 
source to compare it with the models, and therefore predict 



accurately the BH mass and accretion rate for an individual 
object. 



5.3. Do all ULXs in our sample operate in the "same" way? 

Obviously, NGC 55 ULX does not operate in the same way 
as the other objects in the sample. Its large variability am- 
plitude can be explained neither with the case A nor the 
case B possibilities. This is not surprising, given the "dip- 
ping" episodes that have been observed in its light curve, 
and not in other ULXs (or AGN). Due to the presence of 
these "dips" , the source light curves have a large variability 
amplitude for its luminosity. 

IHeil et ail (|2009T ) noticed that the strength of the in- 
trinsic variability in ULXs like NGC 4559 X-l, NGC 5204 
X-l, and Ho II X-l is substantially lower than the observed 
variability amplitude seen in other sources. These objects 
are also included in our sample, but we found that their 
excess variance (i.e. their X-ray variability amplitude) is 
consistent with that observed in the other ULXs. It could 
be due to the technical difficulty of the estimation of the 
PSD fit for faint sources like ULXs. On the other hand, 
we found that cr^xs ^ s n0 ^ constrained in two other ULXs, 
namely NGC 4945 X-2 and NGC 2403 X-l. Given the large 
error on the cr^ xs measurement for these two objects, it 
is not clear whether their intrinsic variability amplitude is 
indeed "significantly smaller" than the amplitude of the 
other sources. In any case though, if the BH mass and ac- 
cretion rate of the objects in the sample are in the range 
of (2.5 - 30) x 10 3 M Q and - 1 - 80%, respectively, we do 
expect their intrinsic excess variance to differ by a factor 
of ~ 100 (sec Fig.[5]). In other words, significant differences 
among the variability amplitude of ULXs do not necessar- 
ily imply that there exist fundamental differences in their 
emission mechanism; they could be due to differences in 
their BH mass and/or accretion rate. 



5.4. The "variability-luminosity" anti-correlation in ULXs 

We found evidence for an anti-correlation between the vari- 
ability amplitude and the X-ray luminosity of the ULXs we 
studied. We estimated the significance of this correlation 
at the 2.5a level. To the best of our knowledge, this is the 
first time that a "variability - luminosity" anti-correlation 
has been detected in ULXs. A si milar anti-correlation was 
also pointed out by IHeil et aP (|2009l ) (see right panel in 
their Fig. 4). They parametrised the variability amplitude 
by means of the power-spectrum amplitude at a given fre- 
quency, which is significantly more difficult to measure 
than the excess variance. As a result, they could measure 
only upper limits on the variability amplitude of many 
sources. Consequently, when excluding the NGC 55 ULX 
data and the points with upper limits, the significance of 
the "variability-luminosity" anti-correlation was substan- 
tially decreased. 

The similarity between the slope of the "variability- 
luminosity" relationship for AGN and ULXs argues in 
favour of the reliability of this anti-correlation in ULXs. 
However, the amplitude of these two relationships is signif- 
icantly different. Under the case B scenario we have dis- 
cussed above, this anti-correlation could be explained if 
the high luminosity objects in the sample have larger BH 
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masses and higher accretion rate. In this case they should 
also have smaller variability amplitudes, as observed. 

However, there are reasons why the correlation shown 
in Fig. [5] may be misleading. Firstly, the luminosity of the 
sources in the sample covers a rather limited range of values 
between (2.5 — 30) x 10 39 erg s -1 . In this case, even a small 
error in the luminosity estimation (due to an inaccurate dis- 
tance and/or flux measurement for example) may shift the 
position of data points in this plot, and hence affect our re- 
sults. We performed a numerical experiment to investigate 
this effect. We used the present (c^ xg — L 2 -iokcv) data set 
(excluding NGC 55 ULX) to create 100 new sets. In each 
run, we randomly decreased and/or increased the luminos- 
ity of all points by a factor of 2. We then fitted the new 
data set (in the log-log space) using ASURV (to take into 
account the upper limits on the data points with negative 
excess variance measurements, exactly as we did with the 
real data points) and recorded the best-fit slope value. In all 
cases, the best-ft slope was different from zero at the (2.2- 
2.5)cr level. The results from this experiment indicate that, 
most probably, "luminosity-induced" uncertainties cannot 
seriously affect the observed anti-correlation. 

However, the most "serious" reason against the relia- 
bility of this anti-correlation is the small size of our sam- 
ple. A larger number of ULXs would be necessary to con- 
firm it. Furthermore, a conclusive test of the "variability- 
amplitude" anti-correlation in ULXs will be possible when 
future observations of ULXs will allow us to determine their 
average excess variance and average X-ray luminosity, and 
examine if they are anti-correlated or not. If this correla- 
tion is confirmed it would indicate that ULXs show a fun- 
damental physical link between BH mass (i.e. luminosity) 
and accretion rate. If this correlation is not confirmed, then 
what we observe in this work can only be a coincidence due 
to incompleteness of the sample. 

6. Conclusions 

The main result of this work is that the variability am- 
plitude of ULXs is significantly smaller than the ampli- 
tude predicted by a simple extrapolation to low luminosities 
of the well established "variability-amplitude" relationship 
for the nearby bright AGN. This discrepancy can be consis- 
tent with the hypothesis that most ULXs operate like AGN, 
but only if: (i) They host an IMBH of - (2.5 - 30) x 10 3 
M Q , (ii) their accretion rate is 1 — 80% of the Eddington 
limit, and (iii) their PSDs have the band-limited noise shape 
shown by GBHs in their "low-hard" and "very-high" state. 

We have also found evidence for an anti-correlation be- 
tween the normalised excess variance and the luminosity for 
ULXs. The slope is consistent with that found in AGN but 
with an offset in luminosity of around four orders of magni- 
tudes. A larger sample of ULXs, and the determination of 
their average X-ray luminosity and variability amplitude, 
is necessary in order to confirm its significance. 
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